Radiative effects in scattering of polarized 
leptons by polarized nucleons and light nuclei 

I. Akushevich a )* A.Ilyichev 6 ), N.Shumeiko & ) 

( a ) North Carolina Central University, Durham, NC 27707, USA 

and 

Jefferson Lab, Newport News, VA 23606, USA 
( b ) National Center of Particle and High Energy Physics, 220040 Minsk, Belarus 

Abstract 

Recent developments in the field of radiative effects in polarized lepton-nuclear scat- 
tering are reviewed. The processes of inclusive, semi-inclusive, diffractive and elastic 
scattering are considered. The explicit formulae obtained within the covariant approach 
are discussed. FORTRAN codes POLRAD, RADGEN, HAPRAD, DIFFRAD and MAS- 
CARAD created on the basis of the formulae are briefly described. Applications for data 
analysis of the current experiments on lepton-nuclear scattering at CERN, DESY, SLAC 
and TJNAF are illustrated by numerical results. 
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Introduction 

Starting with the late nineteen sixties the main source of information about internal 
nature of nucleon is lepton nuclei scattering. Nowadays experimental design of modern 
detectors opens possibility for large set of experimental data during short period of time. 
As a result the calculation of background radiative effects or the procedure of radiative 
corrections (RC) is a quite important thing. 

While considering lepton nucleon scattering it should be noted that the main contri- 
bution to the observables gives subprocess with one photon exchange graphs that is called 
Born process. However all sets of experimental data consist of not only such contribution, 
but of other processes appearing during the measurement which are called background 
ones because of their small contributions to the observables. Therefore the main task 
of data processing in experiments on lepton nucleon scattering is the extraction of the 
one-photon exchange contribution. It leads to the necessity to cancel the background 
events from all data sets. When it cannot be made by experimental methods the corre- 
spond contributions can be calculated in the frame of some theory. Usually a such kind 
of contribution depends on the quantities (the structure function, asymmetries), that are 
observables in the given experiment. As a result we come to the necessity of the itera- 
tion procedure application for the extraction of the observables on Born level from the 
full set of experimental data. This procedure is usually called as the procedure of RC of 
experimental data. 

There are two basic methods of calculation of model independent QED radiative cor- 
rection. First one is connected with introducing of an artificial parameter separating the 
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momentum phase space into soft and hard parts. One can find a classical review intro- 
ducing this formalism in ref. fl],[2|]. However the presence of the artificial parameter is 
a disadvantage of this method. For the soft photon part the calculation is performed in 
such approximation when the photon energy is considered to be small with respect to all 
momenta and masses in the problem. So this parameter should be chosen as small as 
possible to reduce the region evaluated approximately. From the other side it cannot be 
chosen too small because of possible numerical instabilities in calculating of hard-photon 
emission. Calculation within the approach of Mo and Tsai was performed only for the 
case of unpolarized deep inelastic scattering. In the end of seventies Bardin and Shumeiko 
developed an approach || of extraction and cancellation of infrared divergence without 
introducing this artificial parameter. Then this approach was applied to the calculation of 
electromagnetic correction to the lepton current for deep inelastic scattering (DIS) [|| 12 
as well as for some approximation and taking into account the contribution of net soft 
photon emission. The exact expressions for the lowest order QED RC to polarized DIS of 
polarized leptons by polarized light nuclei including targets of spin 1 such as deuterium 
are presented in lljj. Basing on these formulae FORTRAN code POLRAD Q] are con- 
structed. Besides, some others tasks such as Monte Carlo generator RADGEN |l5| was 
made, RC both to diffractive vector meson electroproduction [16,17| and to spin-density 
matrix elements in the exclusive vector meson production |l8| by this method were solved. 
Recently Bardin-Shumeiko approach has being applied to experiment on polarized elastic 

The comprehensive analysis of some results 



electron-proton scattering at TJNAF [19 , 2C 
obtained by BLpl and Q was made in 15, 21|. 

The formulae of || as well as [22] allow along with using the parton model to estimate 
the value of electromagnetic correction to the hadronic current in unpolarized and 
polarized |24|,|^] lepton nucleon DIS. 

Other important source of contributions to systematic error of experiments with high 
energy scattering particles together with electromagnetic effects is electroweak ones. When 
the square transferred momentum is high enough the contribution of Z-boson change is 
compared with electromagnetic one especially on the kinematic borders. Therefore it is 
useful to have the calculation of RC within the standard model of electroweak interactions. 
This task is especially actual for experiments at collider [26,27]. This calculation requires 
choosing renormalization scheme and fixation the gauge. Detailed description of on-mass 
scheme within t'Hooft Feynman gauge were reviewed in [28, 

The calculation of electroweak corrections to the lepton legs for DIS could be found in 
|3C-^. Electroweak corrections to the hadronic current in quark-partonic model (QMP) 
can be found in 1 33 — 36 1 and |[7|] for unpolarized and polarized DIS respectively. 

The calculation of RC within the structure function method was suggest by Kuraev and 



Fadin |38|] and applied for DIS |39| . This method was used for the calculation of Compton 



tensor with heavy photon in the case of unpolarized [40] and longitudinal-polarized 
fermions. The calculation of LO and NLO corrections to radiative tail from elastic peak 
in DIS can be found in [42]. 

The next step is the calculation of higher level RC ~ a 2 . Notice that such kind of 
correction to the elastic tail was calculated in (5|,^2|. As to the continuous spectrum the 
obtaining of the formulae not only in exact form (like J13[| ) but and in ultrarelativistic 
approximation is a quite difficult task (first of all because of the problem of the integration 
over the photon emission phase space). Therefore we consider the correction ~ a 2 to the 



continuous spectrum in the LO approximation. For unpolarized particles such correction 
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in LO and NLO approximation was found in 3S , 53] . 

Another important source of radiative effects in DIS is QCD-correction appearing from 
gluon emission. However QCD corrections, in spite of photon emission is not extracted 
from the structures functions but is directly included into their definition. As a result the 
structure functions start depending not only on a scaling variable x but on the transferred 
momentum squared. 

A lot of articles dedicated to NLO QCD corrections to polarized DIS have been pub- 
lished. The calculations in most of them are performed in the assumption that the quark 
mass is equal to zero (see review |54| and reference there). However there are some works 



when the authors estimate a finite quark mass effect for both unpolarized [46| and po- 
larized structure functions p?|-53|. The main difference of these two approaches is the 
method of tending the quark mass to zero and the procedure of integration of the squared 
matrix element over the phase space of an emitted gluon. In the massless approach a 
fermion mass is equal to zero before the integration while in the massive one this quan- 
tity goes to zero after the integration and survives only in LO terms. 

Notice that QCD and QED radiative corrections (RC) having different origins possess 
some common features on the one-loop level. If our consideration is restricted to so-called 
QCD Compton process then both corrections should be described by the identical sets of 
Feynman graphs. The transition from the strong radiative effects to the electromagnetic 
ones could be performed by the following replacement: 



-a s — * e q aQED- 



(1) 



Due to these facts in [54] the mentioned above Bar din Shumeiko method was applied 
to the investigation of Compton effects within QCD. This calculation is actual because it 
allows us to estimate the influence of the finite quark mass on observables in DIS like it 



was done in [46-51] . 

In the present poster which is basing on ref. |7, 11, 13-pC|, ff4|, |54|] some radiative 
effects in polarized lepton by polarized light nuclei scattering is reviewed. In the first and 
last sections the lowest order corrections are calculated within Bardin Shumeiko approach 
as well as the next order correction for leptonic current are estimated in the leading level 



by the method described in 43]. The calculation of QCD-correction to the hadronic 
current is reviewed with [54|. 



1 Corrections to the lepton current in DIS 
1.1 Inclusive physics 

Let us consider the process of polarized lepton-nucleon DIS 

l(h,0 + N( P , V ) ^l'(k 2 )+X (2) 

where k%, p (£, if) are the momenta (polarization vectors) of the initial particles, h% is the 
final lepton momentum and X is a final hadronic state. 

At the Born or one photon exchange level (fig.[l]a) the cross section of the considered 
process can be presented by the convolution of the leptonic tensor having well known 
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Figure 1: Feynman diagrams contributing to the Born and the model independent QED correc- 
tion of the lowest order. The letters denote the four-momenta and polarizations of corresponding 
particles. 



structure, with the hadronic one Wav The latter can be expanded into Lorentz covariants 
whose coefficients define the structure function which are to be measured. 
As a result the double differential Born cross section has the form [] 

4vra 2 SS X 



0o 



A s Q 4 



j(Q 2 - 2m 2 )9i + (SX - M 2 Q 



2^2\ ^2 



2M 2 



+mMP L {2(Q 2 ^ - qn fc 2 £)^| + (S x k 2 £ - 2 & Q 2 )qvj^) 
+(Q 2 ~ 2m 2 )(Q 2 - 3( OT ) 2 )|| + (SX - M 2 Q 2 )(Q 2 - 3(qn) 2 )^ 

-- ( Q 2 + 4m 2 + 12 nh nk 2 )Q 7 --(Xnk l + S V^QV^j , 

where the kinematic invariants defined in the standard way: 

S = 2k lP , X = 2k 2 p = (1 - y)S, Q 2 = -q 2 = -{k l - k 2 ) 2 = xyS, 

S X = S-X, S P = S + X, \ S = S 2 - Am 2 M 2 , 



(3) 



(4) 



and x = Q 2 /S x , y = S x /S are usual scaling variables. The explicit expression for the 
hadronic tensor and the generalized structure functions 9?j can be found in Appendix A.l 
of®. 

Now and later we deal with both longitudinal and transversal polarized targets. At 
the same time the initial lepton will be polarized longitudinally and its polarization vector 
can be presented as the sum of two parts: 



S 



=fci 



2m 



--P = Co + fi- 



Notice, that in ultrarelativistic approximation over the lepton mass a non-zero contri- 
bution to the polarized part of cross section (||) appears from fo as well as the second part 



1 Now and later a = d 2 a/dxdy. 
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Figure 2: The double differential cross section for inclusive lepton-nuclei scattering sketched at 
a certain value of Q 2 as function of v (arbitrary scale). The three basic channels are the elastic 
(I), the quasielastic (II), and the inelastic (III) one. 



of polarization vector £1 gives a non zero contribution to the cross section of DIS within 
m —> by the correction from real photon emission. 

The natural approximation for high energies of the scattering particles is ultrarela- 
tivistic one 

m 2 ,M 2 < S,X,Q 2 . (6) 
Within this approximation the Born cross section has the form 

4 = - W + (F 2 - ^6 2 )(l - y) - P L P N xy(2 - y) 9l ) (7) 

for longitudinal polarized nucleus and 

CT o = q 4 + ~ bl > x V +( p 2 + — &2)(l-y) 

i " (8) 

-2P L P N ^ Xy ^ M { y 9l + 2 92 )) 

for transversal polarized one. 

It is well-known that there are three scattering channels of virtual boson (7, Z) on 
nucleus in dependence on the transfer energy u = E\ — E 2 , where Ei(E 2 ) is the initial 
(scattered) lepton energy: elastic, quasielastic and inelastic. Representative plot of de- 
pendence of the scattering cross section on u and square of transfer momentum Q 2 = —q 2 
is shown on fig. 1 (only for the regions that give sufficient contribution to RC calculation) . 
the peak in the range I (for v = 0, when neglecting nuclear recoil) corresponds to the 
elastic scattering. Then the nucleus remains in the ground state. Range II stands for the 
quasielastic scattering i.e. direct collisions of leptons with nucleons inside nucleus. A wide 
maximum in the energy spectrum originates from the own movement of nucleons. Range 
III of inelastic scattering occurs when the transfer energy is greater then pion threshold. 

On the Born level both v and Q 2 are fixed by the measurements of the scattered 
photon momentum. Hence, the channel of scattering is fixed too. However, on a level of 
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RC radiated real photon momentum is indefinite, and v and Q 2 are arbitrary so each of 
three channels contributes to cross section. Integration over the photon phase space may 
be presented in the plane of v and Q 2 . Adding the virtual photon contribution a v , we 
have for the RC cross section 

a RC = a in + a d + a< 1 + CT v . (9) 

Here cr m , a el , a q are contributions of radiative tails from continuous spectrum, of the 
elastic scattering radiative tail, of the radiative tail from the quasielastic scattering re- 
spectively. Also the contribution of electroweak correction calculated in the quark-parton 
model is contained in a m . Both a and a 2 corrections are taken into account in a m,el ' q,v . 
To separate the contributions we introduce the lower index, f.e. a m = af 1 + a™. 

A modern approach of solving the task on RC calculation assumes exact calculation 
of the lowest order model independent correction. This correction includes the QED pro- 
cesses of radiation of unobserved real photon and one-loop diagrams to lepton line. These 
processes gives the largest contribution, which can be calculated exactly. Uncertainties 
of the model independent RC can come only from fits and models used for the structure 
functions. The calculation of model dependent correction (box-type diagrams, emission 
by hadrons) requires additional assumptions about hadron interaction, so it has additional 
pure theoretical uncertainties, which are hardly controlled. The model-dependent correc- 
tion is much smaller compared to leptonic radiation because it does not include a large 
logarithmic term log(Q 2 /m 2 ). In this and last sections we concentrate on the calculation 
of model independent correction as main contribution to total RC. 

The explicit expression for the lowest-order QED correction (see ng.|]b-e) can be found 
in [13]. Thus, for infrared free sum af 1 (fig.B(a,b)) and a\ (fig.[l](c,d)) we have 



a{ + a f = *S v a + a f. = + §vert + S l vac + <^ a > + af. (10) 

7T 7T 

The finite contribution from real photon emission 

l(h,0 + N(p, rj) - l'(k 2 ) + 7 (fc) + X (11) 
can be presented in a following way: 



af = - a 3 y 



' min * V 



dR 
~~R 



%(R,t) 3i(0,0) 



(Q 2 + Rt) 2 Q 4 

...ax o 

+ £W j dR WTRW 2 ^ R > T) \> 



(12) 



i= 2 b 



where 



R = 2kp, t = kq/kp (13) 



and & is a momentum of emission photon (k 2 = 0). The limits of integration are defined 

as 

_ W 2 -(M + m n ) 2 _ S X ±^ 

■ftmax — -, , > T ma x min — n , ,9 j 

1 + r 2M A (14) 



X Q = 5 2 + 4M 2 Q 2 , W 2 = S x -Q 2 + M 2 , 
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where is a pion mass. The explicit expression for functions Oij(r) are presented 
in 

The quantity Sj^ appears when the infrared divergence is extracted in accordance with 
the Bardin and Shumeiko method || from a m . The virtual photon contribution consists 
of the lepton vertex correction 5 ver t (figl(d)) and the vacuum polarization (ng.Q(e)) by 
leptons 5 l vac and by hadrons 5^ ac p5| . These corrections are given by formulae (20-25) of 
ref.ll. 

In the case of elastic scattering the nucleus remains in the ground state, so for the 
integration variable R we have an additional relation 

R = R d = (S xA - Q 2 )/(l + r A ), (15) 

resulting in 

1^ a 3 v TA 7 a \ * 4 2MlRiT 2 



,el 

A dxAdy A 2 

TA 



(16) 

Here invariants with the index " A" contain the nucleus momentum p A instead of p (p A = 
M\, Ma is nucleus mass). The quantities Ssf are given in Appendix 1 of 14]. 

Quasielastic scattering corresponds to direct collisions of leptons with nucleons inside 



nucleus. Due to self movement of nucleons we have no additional relation like (|15|). As a 
result we have to integrate numerically both over R and r 



01 



a -i I *E5>M / .m^r^afWr). (17) 



i=lj=l 



The quantities 3f can be obtained in the terms of quasielastic structure functions (so- 
called response functions, see |[4|] for explicit result), which have a form of the peak for 
uj = Q 2 /2M. Due to the absence of enough experimental data the fact is normally used 
for construction of the peak type approximation. The factors at response functions are 
estimated at the peak, and subsequent integration of response functions leads to results in 
terms of suppression factors Se,m,em (or of sum rules for electron- nucleus scattering |56fl): 



3 rmos 4 ki 2M 2 Ri~ 2 

Detailed comparison of the results for quasielastic tail obtained within different ap- 
proaches can be found in 

To take into account the effect of radiation of many soft photons a special procedure 
of exponentiation was applied |0. It can be realized by the following substitutions: 



1 ^ l-xy/A J 1 \ l-xy ) 

OL ( OL \ OL 

(1 + -5 v )a -> exp -5 inf (1 + ~(6 V - 5 inf ))a , 

7T \ 7T / 7T 



1) 

(19) 



ol 

where t r = —(l m — !)• 

7T 
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1.2 Iteration procedure of data processing 

In this subsection the procedure of radiative correction to extract the structure function 
g\{x) from measured spin asymmetry A™ is considered. The spin average structure func- 
tions are considered to be constant and gi (x) equals to 0. The measured asymmetry is 
defined as 

A m = ^ + AA 1 ( JS1 ), (20) 

r\ 

where the radiative correction to asymmetry k.A\ can be written in terms of spin-average 
and spin-dependent parts (a u ' p ) of cross sections @ 

AA _ gffWgi) + + of) - <(gi)K n + < + O (21) 

where 5 V = (Jp/oQ = cZ/ a o- The Born and inelastic radiative tail polarized parts of cross 
sections depend on SF gi, and in the last case the dependence is non-trivial. So the 
equation (^) becomes functional one in g\. This functional equation transforms into a 
system considering the extraction of g\ in concrete binning over x in n kinematic points 
Xi (i = 1, ...,n): 

Af i = ^ + AA 1 (g lj ;j = l,...,n). (22) 

Usually the iteration methods are used to solve such a system of equations. The variant 
of iteration formula is ambiguous, but in practice only two types are used. The first and 
most evident one is to take for fe-th step: 

gf) = F 1 (Af i - AA 1 (g{*-V;j = 1, ...,n)). (23) 

Another possibility to obtain the formulae for iteration procedure arises when both 
Born and radiative correction cross section are separated into spin-averaged and spin- 
dependent parts. Then for the measured asymmetry we have 

Am _ 1 al + a\ _ 91 /F X + of/Doff 



Thus we obtain the iteration formulae 



(k) 

9u 



Af i [\ + 



n u I Drr u 



(25) 



where in right-hand side the dependence on g\ is contained only on the level of RC, but 
not on the Born level. 

The above presented formulae could be applied directly for experimental data on 
hydrogen and deuteron. Three-parameter fit is used for fitting of the spin asymmetry 

IEH§ : 



B/1 a -Cx 



e 



A\(x) = A + x B {l 



Af(x) = {er Ax - 1){B C 



x 



(26) 



For proton asymmetry the one- and two parameter fits are also used: 

A\{x) = x B , A\{x) = Ax B . (27) 
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□ — measured asymmetry 
■ - extracted asymmetry 




ure 3: The results of iteration procedure for spin asymmetry along with the fit constructed. 

For the neutron asymmetry extracted from the deuteron and helium-3 data the nuclear 
corrections have to be taken into account: 



A?(x) 



A 



ZX 3 He 
1 



(x) - (1 - f d (x))P p Al(x) 



fd(x)Pn 



Dilution factors for both cases are given by the formulae 



fd{x) 



1 



+ 1 



fd(x) 



1 



2F V 2 IF% + 1 



(28) 



(29) 



respectively. The numbers P p and P n are effective nucleoli polarizations in nucleus. For 
deuteron target P p = P n = 0.5 — 0.750;^, where tod is the D-state probability (~5%), and 
P p =-0.028 and P n =0.86 for helium-3. 

The application of procedure of RC for experiments on DIS with polarized 3 He and 
HERMES lepton beam with energies E = 27.5 GeV on fig.|| are demonstrated. 

Another important question carefully studied in all experiments on lepton-nucleon 
scattering is contribution to systematical error of measurement due to radiative correc- 
tions. There are several sources for this uncertainty: the RC procedure itself, dependence 
on SF models beyond measured region and physical processes neglected in the standard 
RC procedure. Monte Carlo approach for estimations of these effects is presented in |H| 
for the case of 3 He polarized target. 

1.3 Semi- inclusive physics 

We calculate the radiative corrections to data of semi-inclusive polarized experiments 

e(h) + N(p) — > e'(k 2 ) + h( Ph ) + X( Px ) (30) 

when a hadron h is detected in coincidence with the outgoing lepton. In this case the 
cross section depends on five kinematic variables which can be chosen as 



x, y; z, t, 4> h , 



(31) 
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where x and y are usual scaling variables, z and t are defined via hadron momentum 



t=(q-Ph) 2 , z=p h p/pq, q = k 1 -k 2 , (32) 

4>h is an angle between planes (fci,^) an d (q,Ph) m the rest frame (p = (M, 0)). In- 
stead of t— dependence we will also consider the cross section as a function of transversal 
momentum pt defined in (A. 8) of [jll|]. Also the following invariants will be used 

M 2 X =p 2 x = (l- z)S x + M 2 + t, V lt2 = 2k 1>2Ph . (33) 

Notice that there are two tasks for this process. The first task corresponds the de- 
scription of semi-inclusive DIS in the frame of quark-parton model. In this way the cross 
section depends on three variables x, y and z. However for the real situation the inte- 
gration region is limited by the experimental cuts. To deal with them on the angles of 
registered hadron a special procedure was developed and cross section describing such 
situation is defined by five variables: three standard x, y, z and two additional pt, 4>h- 
The relation between these two cross sections can be established as follows: 

d&R f 1 d?k dz . 2 .9 a\n/ 2 -2a ^CT 



/l d k dz 
— — — — dp t d(f> h 6(sm 2 -draax - sin 2 i9)6»(sm 2 <& - sin 2 i? min ) 
2tt 2kn dz 



dxdydz J 2ir 2k$ dz dxdydzdptdcpu 

(34) 

if ^-function are removed. Here $ is the angle between the beam direction (ki) and hadron 
momentum in the lab frame. 

Within quark-parton model for the Born cross section of longitudinal polarized semi- 
inclusive DIS we have 

where 

F U = 2(1/|/ - 1 - fi N x) + y, FP = y-2, m = M 2 /S. (36) 
The quantities 

£+H( Xl *) =^e 2 [/+(x)±/-(x)] J Df (z) (37) 
g 

is the ordinary for QPM combination of the distribution functions fq^ \x) for the quark 
of flower q polarized (anti)parallel to the nucleon polarization, and of the fragmentation 
functions D^(z) of the quark q into the hadron h, e q being the quark charge in units of 
elementary charge. 

As it was mentioned above there are three radiative tails give contribution to 
DIS:inelastic, elastic and quasielastic (see Fig ^). However, in semi-inclusive process 
the transfer energy is above pion threshold, so in RC calculation one have to take into 
account only a vn and a v . The lowest order QED correction was calculated in ref. and 
presented in []l4| as the sum of factorizing and non-factorizing infrared free parts 

_ d 3 a E M at, F 
aEM= dx^yTz = ^ VRao + aR > (38) 

where the first one is the sum of infrared divergence separated from radiative tail and the 
additional virtual photon contribution: 5]/r = 5y + 5 R R . 



11 



The finite part of (pq) has the form 

°R = -^r [dh [ dt 2 {^-[F^+(x,z) + P L P N F^-(x,z)]- 

S I I XtVt (39) 

J{tu- h) [F u rE+{Xj z) + PlPnF p r j:- {Xj z)]h 

xy 

The definition of the rest quantities can be found in 
In order to write five-dimensional cross section it is useful to take the some scalar 
products with ph in the form 

2^1,2 = h,2Ph = a 1 ' 2 + b cos <j> h (40) 
-^i? = kph = R(a k + b k (cos <ph cos <pk + sin cfih sin <pk))- 

where fi = kph/kp and 4>k is the rest frame angle between planes (ki,k 2 ) an d (q,k). We 
will use = a 2 ± a 1 . The explicit expressions for coefficients are given in 

As a result the five-dimensional cross section of unpolarized DIS on the Born level has 
the following dependence on tp^'- 



a ° = dxdvdzdr> 2 d6 h =^ {A + C ° S ^ + C ° S 2<PhAC% ' " ' 



where N = a 2 yS x / ^/Aq. The coefficients A, A c and A cc do not depend on 0^ more and 
they have the form 

A = 2Q 2 H 1 + (SX-M 2 Q 2 )H 2 + (4a 1 a 2 + 2b 2 -MlQ 2 )H 3 

+ (2Xa 1 + 25a 2 - zS x Q 2 )H^ 
A c = 2b{2a + H 3 + S P H 4 ), 

A cc = 2b 2 H 3 . (42) 

The model for structure functions Hi can be constructed on basis of results of the paper 
Mulders and Tangerman [Q. Keeping only the leading twist contribution we have for 
structure functions 

2 2 

« ^ i (43) 

and 

g = g 1 = beM-bPt), (44) 

where b = R 2 /z 2 is a slope parameter and R is a parameter of the model. Alternatively 
power fit is also considered in [11]. 

When integrated over the kinematic variables <f>h and pt the cross section ( |4"l"| ) coincides 
with the unpolarized part of semi-inclusive cross section calculated within quark parton 
model p|. 
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The cross section that takes into account radiative effects can be written as 

cr bs = &oe S ' nf (1 + 5vr + Svac) + <r F . (45) 



Here the corrections 5i n f and 5 vac come from radiation of soft photons [12| and effects of 
vacuum polarization]^. The correction 5vr is infrared free sum of factorized parts of real 
and virtual photon radiation. These quantities are given by the following expressions 



a 



3, n 1 2 X> S'X'-Q 2 p 2 x it 2 



5vr = —[ —l m —2 In h Li2 

vh k\2 2 S' S'X 

A a n (Pi ~ (M + m n ) 2 ) 2 

Oinf = -(im-l)ln -^-^ , (46) 

e fd e Vt _i_ 

°vac — °vac °tiac j 



tiac ' 

where S" = X + Q 2 - V 2 , X' = S - Q 2 - V\, l m = In Q 2 /m 2 and Li 2 is Spence function or 
dilogarithm. 

The contribution of radiative tail has the standard form (see eq.(|l2|)) 



2"7T T-max A *^ Rmax 

t i=l j=l n 



Hi , H 



(Q 2 + i?r) 2 J Q 4 



.(47) 



Here 2M 2 W . min = S x ± y^Q and = (M 2 - (M + m n ) 2 )/{l + r - /i), <5 i= l for j=l 

and 5j=0 otherwise. The explicit formulae for functions 9(t, cpk) can be found in [11]. The 
structure functions have to be calculated for Born kinematics, but TLi is calculated in 
terms of tilde variables introducing by 

Q 2 = ( q - k) 2 = Q 2 + Rt, 

W 2 = {p + q-k) 2 =W 2 - R{1 + t), 

t = (q-k-p h ) 2 =t + R(T-fi), 

M 2 X = p 2 x = Ml + R(l + r-fi). (48) 



1.4 Higher order effects 

The dominant contribution to RC can be obtained by the method of leading logarithms. 
Being based on the fermion mass factorization this approximation allows us to calculate 
corrections ~ {a\og{Q 2 / m 2 )) n . The below formulae were obtained in [Q by the method 
which is presented in |]43| . 

The contribution of inelastic tail a % 2 together with loop effects a 2 

—in i „v tin _ i „in , Jin i —in , in , in , in , in , in , in fAri\ 

a 2 + a 2 = d 2 a + a Vkl + a Vk2 + o klkl + o k2k2 + a klk2 + a lkl + a lk2 +a fkl +<? fk2 . (49) 
The factorized part is 

82 = iJ( 3< W) + 2l mSs P 5(Q 2 ) + + 4Li 2 (l - z 2 ) + 12Li 2 (l - Zl ) - ^vr 2 ]), (50) 

2 There are explicit formulae for leptonic contribution to vacuum polarization effect (see ]l3[ for example) 
and parameterization of hadronic one [p5| . 
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l m = logQ 2 /™ 2 , S v = 21og((l - Zl )(l - z 2 )) + 3, and 5{Q 2 ) = S l vac + <5£ ac . 

The contribution from the vacuum polarization if it coincides with real photon radia- 
tion has the form 



a Vki,Vk 2 



a 2 f dz 

^2 lm J — 

21,2 



(1 + z 2 )5(t X:S )a s , p - 25(Q 2 )a 



(51) 



where t x = zQ 2 and t s = Q 2 /z. 

Next three terms correspond to the cases when two radiated photons are collinear to 
an incident electron (cr klkl ), an outgoing electron (a k2k2 ) 



- —l 2 f A, 

a kik 1 ,k 2 k 2 ~ 87T 2 m / 

zi, 2 



— (21og[(l - z)(l - z 1>2 (z))] - lnz + 3)x 



x ((1 + z 2 )o kx2 - 2a ) + ((1 + z) lnz - 2(1 - z))a fel2 



(52) 



or when one photon is radiated in incident electron direction and the other in the outgoing 
electron direction {a klk2 ): 



in 
(T k 1 k 2 



Of 

47 2 ' 



2 fjz_ I _ 
m J \ -z> J 1 

z i z 2 {z') 

-2(1 + z 2 )a k2 + Aa 



dz 



(1 + ^)(i + z ' 2 )a klk2 - 2(1 + z'> fcl 



I2\ 



(53) 



Here 



Vk lk2 = yMz'S, X/z, z'Q 2 /z)/(z(zz' - 1 + y)), 
zi(z) = (1 - y)/(z - xy),z 2 (z) = (l-y + zxy)/z. 



(54) 



There are two channels (singlet and non-singlet) of the fermion pair production that 
give a contribution to a 2 of order RC. The singlet channel 



l-4m ; iW/S 



in _ ^ 72 
a lk 1 ,lk 2 — ^~^ L 



a 

8^ 2 " m 



dz(2(l + z) log z + 1 - z + -(1 - z 3 )z)a kl 2 (55) 



21,2 



corresponds to the case when the incident and the outgoing lepton as well as the leptons of 
the unregistered pair belong to different leptonic lines connected by an additional virtual 
photon. 

The rest non-singlet part 



l-4m f M/S 



21,2 



(1-z) 



(56) 



arises from the two-lepton decay of an additional virtual photon. 
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The main contribution to the second order elastic and quasielastic radiative tail arises 
when the additional radiated photon is collinear to the incident or outgoing fermion line: 



°f = 7^ l m$s P o-f + af lt + af 2t , (57) 



where 



el « 



(1 + z A )af - 2af 
< a t = J dz"- Lh* (58) 



The quantities cr£' are obtained in the terms of approximate elastic radiative tail af = 
af(x,y,S): 

< = yo?(xyz/(z + y-l),(z + y- l)/z, zS)/(z -l + y), 

(59) 

< = yaf(xy/(z + y -l),( z + y - l)/ z , S)/z(z -l + y). 

The numerical analysis shown that correction ~ a 2 to the considered asymmetry is 
about some per cent from the full RC and has a maximum near kinematic borders of the 
considered experiments. 

However it is necessary to remind that correction ~ a 2 has been calculated within 
LO approximation while even ultrarelativistic approximation for radiative tails from 
elastic and quasielastic peaks leads to high corrections. Besides, the approximation 
~ a 2 log 2 (Q 2 /m 2 ) gives rather good results for the factorizad part of correction but it 
produce poorer ones for the remaining components which make the dominant contribu- 
tion to the asymmetry in polarized particle scattering. Hence, a more accurate calculation 
is required for more detail analysis of ~ a 2 correction. 



1.5 FORTRAN codes POLRAD 2.0 and HAPRAD 

Using the results presented in the previous subsections two FORTRAN code POLRAD 
2.0 [|l4j] and HAPRAD have been created. Here some results obtained by these codes 
are presented. 

The program POLRAD 2.0 consists of two patches. The first one which itself is called 
POLRAD describes RC to inclusive DIS. In order to estimate some radiative effects by 
this patch let us separate cross section into spin-averaged and spin-dependent parts 

0-0,1,2 = 00,1,2 + p N°a,i,2- (6°) 
and consider the following quantities: 

u,p . u,p u,p 1 u,p 1 u,p 

x u,p _ °0 + a \ cu,p _ °j + a \ + cr 2 / R1 \ 

°1 — u,p j °2 — u,p ) V 01 / 

which are characterized the ratio of the total unpolarized and polarized part of the cross 
section to the Born one (with and without second order effect). The dependence of these 
quantities on scaling variable in HERMES kinematics is presented in fig. ||[ 

The second patch of POLRAD 2.0 which is called as SIRAD allows us to calculate RC 
to three-dimension cross section da / dxdydz of semi- inclusive DIS. Fig. || shows the ratio 
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a) b) 

Figure 4: x— dependence of the ratio of the total cross section to the Born on with (dashed 
curves) and without (solid curves) the second order effects for electron-proton DIS with HER- 
MES kinematics at the initial lepton energy of E e = 27.5GeV for 1) y = 0.1, 2) y = 0.3, 3) 
y = 0.6, 4) V = 0-8/ Q>)the correction to unpolarized part of cross section; b)the correction to 
polarized part of cross section. 
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Figure 5: The radiative correction factor to semi-inclusive cross sections (r) = cr/cr , upper 
plots) and to polarization asymmetry (5a = (A obs — A born )/ A born , lower plots) versus z. The 



cross sections were integrated over other kinematic variables within HERMES cuts 
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X 


y 


Q 2 

GeV 2 


z 


SIRAD 


HAPRAD 




without cuts 


with cuts 


Qi 


Q'x 


Q2 


Qx 


Si 


Q2 


0.038 


0.677 


1.33 


0.25 


1.029 


1.033 


1.024 


0.982 


1.041 


1.025 


0.985 


0.062 


0.567 


1.82 


0.35 


0.996 


0.989 


0.989 


0.947 


0.989 


0.980 


0.951 


0.092 


0.529 


2.52 


0.45 


0.970 


0.961 


0.961 


0.934 


0.961 


0.956 


0.936 


0.131 


0.499 


3.38 


0.55 


0.945 


0.936 


0.933 


0.912 


0.934 


0.931 


0.906 


0.198 


0.476 


4.88 


0.65 


0.918 


0.902 


0.902 


0.889 


0.897 


0.897 


0.881 



Table 1: The results for RC factors to three dimensional semi-inclusive cross section obtained by 
FORTRAN codes SIRAD and HAPRAD (see text for further explanations). Kinematic points 
are taken from the Table 1 in ref. ffiBJ . 



r] = a obs inborn an( j relative correction to asymmetry 6 a = (A obs — A born ) / A born for K 
and 7r mesonproduction. 

To estimate the RC to five-dimensional cross section d 5 a/dxdydzdpfd6h the special 
FORTRAN code HAPRAD was developed. Now we show that numerical results for RC 
to da I dxdydz reproduced by these two codes which coincide with a good accuracy. It 
can be seen in the Table [l] where we represent RC to the cross section as it follows from 
the runs of the codes POLRAD 2.0 and HAPRAD. Since HAPRAD allows us to take 
into account the kinematic cuts and to use different models for p^-slope, three fits for 
p\ -distribution and the cases with and without experimental cuts were considered. The 
first fit for pi -slope is defined in (44), while the second and third ones are our fits of 
experimental data [64| using exponential (Q[ = Q\ at b = a/z) and power functional 
forms 

1 



Q = Q 2 



c+dz 



(62) 



a + bz + pi _ 

As the kinematic cuts on fa and pt of the measured hadron we took HERMES geomet- 
rical cuts We can conclude from this analysis that neither important differences 
between SIRAD and HAPRAD results, if exponential model for ^-distribution is used, 
nor dependence on slope parameter model and applying of geometrical cuts are found. 
However RC takes a negative shift in the model based on power functional form (|62|), As 
it is shown in [|llj RC depends on steepness of pi distribution. It is a reason why models 



like 6(pf) (QPM, POLRAD 2.0) and (gg) give larger RC. Within practical RC procedure 
in concrete measurement of da /dxdydz the model can be fixed only if the information 
about p^-distribution is additionally considered. 



1.6 Monte Carlo generator RADGEN 1.0 

In this subsection we present a Monte Carlo generator RADGEN 1.0 ]l5| for the events 
with a possible radiation of a real photon in DIS on polarized and unpolarized targets. 
The events are generated in accordance with their contribution to the observed total cross 
section given by 

&obs = Cnon-radi^) + °m(A) + a q + a d . (63) 
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The first term (T n on—rad(A) contains not only the contribution from the Born process but 
the contributions from loop corrections (a v ) and from multiple soft photon production with 
a total energy not exceeding a cut-off parameter A. In handling the radiative corrections 
two approaches have been used, the one developed by Mo and Tsai ||],||] and the other 
given by Bardin and Shumeiko ||. 

Kinematics of the Born (or one photon exchange) process is completely defined by 
the scattering angle 9 and the energy E' of the scattered lepton, two variables which are 
usually measured. All other inclusive kinematic variables can be expressed in the terms 
of 9 and E' . Having in mind the laboratory frame, i.e. the frame with the target nucleus 
at rest the relevant variables are 

9 O 2 
Q 2 = AEE' sin 2 -, u = E - E' x = -f— , 

2 2M,' (64) 

y = ^, W 2 = Q 2 (l/x-l) + M 2 , 

where E is the beam energy. 

The event registered in the detector with certain values of E' and 9 for the scattered 
lepton can either be a non-radiative or a radiative event, i.e. an event containing a 
real hard radiated photon. For the radiative event there are three additional variables 
necessary to fix the kinematics of the radiated real photon apart from E' and 9. A 
possible choice is the photon energy E 7 and the two angles # 7 and 7 , where 6* 7 is the 
angle between the real and the virtual photon momenta k and q = k\ — &?2 and (j)-y is the 
angle between the planes defined by the momenta (k\, fo) and (k, q). For the events 
with the radiation of a real photon the kinematic variables describing the virtual photon 
and which is used to generate the hadronic final state is derived from eq.(^4|) because the 
substitution q — > q — k has to be made in their definition. The variables obtained after 
this substitution will be referred to as 'true' ones: 



Wt rue = W 2 - 2E 1 (y + M - \jv 2 + Q 2 cos 7 ), v true = v - E y , 

Qt 



Qtrue = Q l + 1E^{v- Jl/ 2 + Q 2 COS0 7 ), X t 



'true 



2Mu tr 



(65) 



For non-radiative events the true kinematics exactly coincide with (|64|), 

In case of the elastic and quasielastic radiative processes the mass squared of the 
hadronic final state is fixed imposing additional constraints on the true kinematic vari- 
ables. Indeed the following ranges are allowed: 

x < xtme < 1 for a in 

Xtrue = 1 f° r &q Qrnin — Qtrue — Qmax (66) 



xtrue = M A /M for a, 



el 



where 



r> 2 _ n2 2(l-x r )(l±v/lT^)+7 2 , , 
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and 

7 2 = 4M 2 x 2 /Q 2 , X r = x/xtrue- 

The Monte-Carlo procedure for the generation of events with a possible photon radi- 
ation is the following: 

The procedure starts off with the generation of the kinematics of the scattered lepton 
and the calculation of an event weight from these kinematics. Then the appropriate 
scattering channel (non-radiative; elastic, quasielastic or inelastic radiative tail) has to be 
chosen according to their contribution to the total observed cross section (see (|6^)). If a 
radiative channel is selected the radiated photon has to be generated and the values of the 
kinematic variables have to be re-computed to obtain the true values. For each event the 
weight has to be recalculated. The new weight is defined as the ratio of the radiatively 
corrected and the Born cross section. After this recalculation the weighted sum of all 
events (generated originally in accordance with the Born cross section) gives the observed 
cross section. 

The recalculation of the weight requires the knowledge of the cross section integrated 
over the photon momentum. This is done differently in two theoretical approaches. In 
the Mo-Tsai approach an additional parameter A dividing the integration region into a 
soft and a hard photonic part is introduced. There is no such parameter in the Bardin- 
Shumeiko approach. However, in both approaches a minimal photon energy EZ- is 
adapted in generating a radiated photon. If the photon energy is above this value the 
kinematics of the photon is calculated and the event becomes a radiative one. The actual 
value of depends on the aim of the photon generator. In general photons should be 
detected in the calorimeter. So the energy threshold of the calorimeter sets the value for 
E 1 ■ . 

mm 

The unpolarized generator is based on the FORTRAN code FERRAD35 |6(| . The code 
calculates the radiative correction to deep inelastic scattering of unpolarized particles in 
accordance with the analytical formulae given by Mo and Tsai @,g. 

The result for the lowest order radiative correction of the cross section is 

o- = 5r(A)(1 + 5 vert + 5 vac + <5 sm )<7i 7 + (T e l + Vq + fin (A), (69) 

where <ti 7 is the one photon exchange Born cross section and S vac , 5 ver t, and 5 sm are 
corrections due to vacuum polarization by electron and muon pairs, vertex corrections 
and residuum of the cancellation of infrared divergent terms independent of A (see [|nj 
for detail). 

The cross sections a e i, cr q , and Oi n are the contributions from radiative processes 
(i.e. including real photon radiation) for elastic, quasielastic and deep inelastic scattering 
(see Hflim for detail). 



An artificial parameter A had to be introduced to divide the integration region over the 
photon energy into two parts, the soft and the hard energy regions. The hard energy region 
tri n (A) can be calculated without any approximations. The soft photon part is calculated 
for photon energies approaching zero. After cancellation of the infrared divergences and 
a resummation of soft multiphoton effects the correction factor is given by 



5 R (A) = exp 



^(inf + mfm^-, 




(70) 
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The polarized generator is constructed utilizing the FORTRAN code POLRAD 2.0 



1 13, 14 1 that was presented above. The cross section formulae (|9r]l q) can be rewritten into 
a form similar to eq. (|69|) : 



° = S R( E min)( 1 + S vert + <W + S sm )a ll + (7 add (E^ in ) + cr el + a g + <r in (£^ n ). (71) 

In the ultrarelativistic approximation (Q 2 3> w?) the corrections 5 ver t, S vac , and 5 sm 
take the same form as in eq.(|69|) 0. Since this code is based on the method of covari- 
ant cancellation of infrared divergences developed by Bardin and Shumeiko 11] the final 
formula for the cross section free of infrared divergences does not include any artificial 
parameter like A. However, ET- has to be introduced if the code is used considering the 
generation of real radiated photons. As a result the formula (|l2|) has to be rewritten in 
the following way 

Traax 4 Rinax _ 2 

o m = -o?y I drX)E%W / dR (Q2 + RT) 2 %{R > T) > (?2) 
. j=i 7=1 R / 

while the elastic and quasielastic radiative tails were defined above by eq.(|l6|) and ( |l8| ) 
respectively. The lower integration limits is determined by R m in = ^ME 1 ^^. 

Apart from the usual ultrarelativistic approximation another one was made when the 
eqs. ( p9| ) and ( |3l| ) were obtained, namely the photon energy was considered to be small 
E 1 <C E, E' in the region E~* < W min . In POLRAD, the term a add {EP min ) is added to take 
the difference of this approximation to the exact formula into account: 



'max 4 1 

t r i=l n 



^(i?,r) 9,(0,0) 



(Q 2 + i?r) 2 g 4 



2Af£ 7 - 

mm -,7-2 



~, Rj z 




[73) 



Therefore, the results obtained with the help of POLRAD are independent of the thresh- 
old parameter E^ min by construction. 

To estimate the value of radiative effects numerically Monte-Carlo event samples of 
100k events each have been generated for a 3 He target and kinematic cuts relevant for 



the HERMES experiment |7| have been applied (Q 2 >1 GeV 2 , W 2 >4 GeV 2 , y <0.85, 
0.037< 9 <0.l4rad). The distributions of the radiation angles 6? 7 and 7 as defined in 
the so-called Tsai-system in which the z-axis is along the direction of the virtual photon 
and the y-axis is normal to the scattering plane (see ref. Jl],|2)) and of the energy E~ are 
displayed in fig.||. These distributions are generated for a certain point in the kinematic 
plane of the scattered lepton, i.e. for x = 0.1 and y = 0.8. The two peaks seen in the 
# 7 distribution of fig|3| correspond to the s- and p-peak emerging from collinear radiation 
along the direction of the incoming and outgoing lepton, respectively. 



J It should be noted that contributions from r-leptons and quark loops were included into 5 V 
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Figure 6: The distribution of the radiation angles 9 1 a) and 7 b) and of the energy c) of the 
radiated photon for x — 0.1 and y = 0.8. The two-dimensional distribution d) shows 9 1 vs </> 7 . 



1.7 Electroweak correction 

It is natural, that data processing of modern experiments on DIS requires correct account 
of RC. The above presented expressions for RC to polarized DIS within QED is successfully 
used for fixed target experiments. 

However when the polarization DIS experiments at collider will be possible p7|| , we 
cannot restrict our consideration to 7-exchange graphs only, because in this case the 
squared transfer momentum Q 2 is so high that weak effects begin to play an essential 
role in the total cross section and spin asymmetries. At the same time the ratio m 2 /Q 2 
(where m is the mass of a scattering lepton) becomes so small that it could be restricted 
to non-vanishing terms for m — > 0. The such calculation was already done by us [37|, but 
there we used only the naive parton model. 

In this subsection the electroweak correction to the lepton currents within QCD- 
improved parton distribution. 

Generally, the double differential cross section lepton-nucleon scattering in the frame 
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electroweak theory has a form 

! (74) 



4ira 2 S x S 

cr Q : 



XsQ 4 



LJlW^ip, q) + \{Lf v + L$)Wtf(p, q) X L Z ^W^ z (p, q) X 2 



where x = Q 2 /(Q 2 + M§) (M z is the Z-b oson mass), while the explicit expression for 
leptonic L l ^ v and hadronic Wff v tensor = 7, Z) can be found in [32|. 



After tensor contractions cross section ( |74|) can be written in more simple form 

4-Ka 2 y 



^ i=i 



Y^OfFi (75) 



through linear combination eight generalized structure functions: 

ft = RlF? + xRfFf + X 2 R^F? (i = 1, 2, 6 - 8), 

(76) 

Ti = R\F? + X R\ Z Ff + X 2 R Z A F? (i = 3 - 5). 
Here the quadratic combinations of the electroweak coupling constants are defined as: 

m n = (v m v n + a m a n ) - P L (v m a n + v n a m ), 

(77) 

R mn = ( v m a n + _ p^m^n + 

where 

= 1, w z = (-1 + 44)/4s wCu „ , . 

aT = 0, a z = - 1/4S,,, c,,,, 1 ] 

c w and are cosin and sine of Weinberg's angle respectively and Pl is degree of lep- 
ton polarization. The quantities Of depend only on the target polarization vector and 
kinematic invariants: 

Of = Q 2 , 0§ = rjqQ 2 S p /2M 3 , 

Of = (SX - M 2 Q 2 )/2M 2 , Of = -{Xrjkx + Sr/k 2 )/2M, 

Of = Q 2 S P /AM 2 , Of = m(SX - M 2 Q 2 )/2M 3 , 

Of = -Q 2 v(h + k 2 )/M, Of = -r]qQ 2 /M. 

Model independent part of lowest-order electroweak correction which includes contri- 
butions both from a real photon emission (<tr), and from the additional virtual particles 
(ay) can be presented as a sum of infrared free terms: 

<?RC = OV + a R = —fivRCQ + °V + a R + ( 80 ) 
7T 



The factor 



A n Q2 1M {W 2 -{M + m^) 2 ) 2 3, Q 2 

5vR = (1 ° g ^ " 1)ln (X + Qi)(S-Q 2 ) + 2 lQg ^~ 2 



1 2 X + Q 2 . SX-Q 2 M 2 vr 2 
-o lo S -s^Q2 +Ll2 {X + Q 2 )(S-Q 2 )~~6 M) 
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appears in front of the Born cross section after cancellation of infrared divergence by 
summing of an infrared part separated from an and that part of virtual contribution 
which arises from the lepton vertex graphs including an additional virtual photon. 

The contribution from electroweak rest loop correction can be written in terms of the 
Born cross section with the following replacement: 

a r v = a B - 5RV n A ) , (82) 

and SR^JX defined in ||. 

The infrared free part of the cross section of the radiative process can be presented 
as a sum of two terms <rjj and gr. The first term includes unpolarized part of the cross 
section as well as that part of cross section which appears from leading term £o of lepton 
polarization vector (||). It looks like eq. 

^(0,0) " 

Q A _ 

(83) 

where the integration variables and their limits are defined by formulae flHSynp . The sec- 
ond term containing the contribution which is proportional to the rest part of polarization 
vector £i in the limit m — > can be written through the Born cross section: 

^max 

ay f Ran B 

aR = ^s J W^R) apl 

The upper integration limit defines as R s max = S(W 2 — (M +m 7T ) 2 ) / (S — Q 2 ), and aj^ is the 
proportional Pl part of the Born cross section with the following replacement kinematic 
variable: S -> S - R, Q 2 -> Q 2 {1 - R/S) and hrj -> hr](l - R/S). 

The double differential cross section as a function of the polarization characteristics 
of the scattering particles can be presented as the sum of four terms: 

a = a u + P L <jZ + PncjV + P N P L <jto, (85) 

the first of them is an unpolarized cross section and three others characterize the polar- 
ized contributions independent on polarization degrees. There are no problems with the 
luminosity measurement in the current collider experiments, so apart from the usual mea- 
surement of polarized asymmetries the absolute measurement of the cross sections with 
different polarization configurations of beam and target will be probably possible in future 
polarization experiments at collider. Besides, now the new methods of data processing, 
when experimental information of spin observables is extracted directly from the polarized 
part of the cross section p8|, |6^| are actively developed. In it is shown how to sepa- 
rate completely unpolarized and polarized cross sections from a sample of experimental 
data using a special likelihood procedure and a binnig on polarization degrees. All above 
mentioned as well as the fact that RC to asymmetry is always constructed from RC to 
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a 3 y / dr^Wir) 



Rmax 



i=l 

Rmax 



dR 
~~R 



Ti(R,T) 



+E%( r ) / dR 



i=2 



\Q 2 + Rt) 2 



(Q 2 + Rrf 



Fi(R,r] 



(84) 




Figure 7: Radiative correction to Act defined in (ffiW with (dashed curves) and without cut (full 
curves). The down indexes \\ and _L correspond to longitudinally and transversely polarized 
proton beam respectively. 
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parts of the cross section allows to restrict our consideration to numerical studying of RC 
to all of the cross sections in the equation ( |S5|) and their combinations. 

Radiative corrections to these cross sections defined as a ratio of the cross section 
including one-loop RC only to the Born one 

$ a = % = (a = u,C,V,^) (86) 

is presented on fig. [7| as a function of scaling variables x and y. The kinematic region 
corresponds to the present unpolarization collider experiment at HERA J26[ |. 

In experiments at collider a detection of a hard photon in calorimeter is used to reduce 
radiative effects. The dashed lines on fig. ^ demonstrate the influence of an experimental 
cut on the RC. One of the simplest variant of the cut when events having a radiative 
photon energy EL > lOGeV are rejected from analysis, is considered only. 



2 Correction to the hadronic current in DIS 
2.1 Tools of calculation and explicit results 

On the language of the hadronic tensor the lowest-order one-loop RC to the hadronic 
current consists of two parts whose contributions are calculated in a different way: 
^}w lo ° v = W R U + W^ u . The first one appears from a gluon emission and requires the 
integration over its phase space, while the second part comes from a gluon exchange 
graph and reads 



3 7T Z - 

q 



^ + ^)c-«-^ + 5'.- J+ T 



wS! + W* MM , (87) 



where is a contribution of g-quark to the hadronic tensor on the Born level, and 
Wf^' IM is the quark anomalous magnetic moment. The pole term which corresponds to 
the infrared divergence is contained in V IR . The arbitrary parameter \x has a dimension 
of a mass. 

Both of these contributions include the infrared divergences, which have to be careful 
considered in order to be canceled. Like QED we use the identity 

W R = W R - W IR + W IR = W F + W IR 

Here W F U is finite for k — > 0, and Wl R is the infrared divergent part of W R V . Using the 
dimensional regularization scheme the latter can be given in the form 



4a y 



3 - , 



2 (^ H + ln^a)(i g -l)+l^ + ^ 

w 0q 



- 7 -i + 3 --^ 

2 l v 4 Lq 4 v + 4 3 



where l q = \og(Q 2 / m 2 q ) , l v = log((l - x)/x). The sum of WfJ} and W^ u 

q 

3 n l—i u q yv ~ r w 
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is infrared free. 

In order to extract some information about QCD contribution to the polarized struc- 
ture functions, the integration in W^ v over the gluon phase space should be performed 
without any assumptions about the polarization vector rj. So the technique of tensor in- 
tegration have to be applied in this case. Since the result of the analytical integration has 
the same tensor structure as the usual hadronic tensor in polarized DIS, the coefficients 
in front of the corresponding tensor structures (like g^vi P(j,Pv ■■■ 

) can be interpreted as 
one-loop QCD contributions to the corresponding structure functions. 
Thus the QCD-improved structure functions -Fi,2,l and gi read 

Fi(x, Q 2 ) = ^[F2(x, Q 2 ) - F L (x, Q 2 )], F 2 (x, Q 2 ) = e 2 q f q (x, Q 2 ), 



Fl(x,Q 2 )-- 
and g2 looks like 

92(x,Q 2 ) 



-^- x J2 e q [ dz fq( x / z )' 9i{x,Q 2 ) = \^e 2 q Af q (x,Q 2 ), 

°" i i i 



;i-2Z ff -ln(l-x))A/ ff (z)+ / dz (41. 



-41ogz(l - z) - 12 



(1-z) 



(90) 



where the <5 2 -dependent unpolarized and polarized parton distributions are defined as 



2a, 



2a, 



dz\,l + z l 



f q (x, Q 2 ) = (1 + -^5 q )f q (x) + -±J -[(- 
7 1 



2 1 



_ + 3* + 4)/,H 

z z 

2a 



2a, 



l q + log 



1 



[l q -\ogz{l-z)) 



t ) /«(*) 



AUx, Q 2 ) = (1 + -±5 q )Af q (x) + -± / - ( 



3tt 



3tt 
2 



r.l + z' 



1 



lq + log 



2 
Z 



(l q -logz(l-z)) 



1 



4 ) ^) 



(91) 



and <5 g can be found in (|88|). 

It has to be noted that our formulae (Hl^iD are in the agreement with ones obtained 
earlier. The unpolarized structure functions coincide with (2.24,2.49) from [46], gi(x,Q 2 ) 
corresponds to the expression (13) in [47]. At last g2(x,Q 2 ) can be considered with the 
sum (18) and (19) from @. 

After integration of the QCD-improved structure functions over the scaling variable x 
it could be seen that QCD RC to the first moment of the unpolarized structure functions 
as well as g 2 have the identical values both for massive and massless quark approaches. 
However the value of QCD-correction to the first moment of g\ depends on the scheme. 
It was calculated: 



gi(x,Q 2 



7T 



dxg1(x), 



(92) 



where C g \ = 1 for massless and C g \ = 5/3 for massive one. 
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2.2 Discussion and Conclusion 



In this section we applied the approach traditionally used in QED and electroweak theory 
for calculation of QCD correction to the DIS structure functions and sum rules. Since the 
quark was considered massive within this approach, it allowed us to estimate the finite 
quark mass effects at the NLO level. LO correction contributes to the DIS structure 
functions but vanishes for the sum rules, so this mass effects are important just for the 
sum rules. 

We found that there is no any additional effect for the first moment of the unpolarized 
structure functions as well as for Q2 ■ However there is some non-zero correction to the first 
moment of g\ and as a result to the Ellis- Jaffe and Bjorken sum rules. The value of the 



correction is in agreement with the results of refs. |47, 5C||, obtained by different methods 



We confirm also the statement of [47| that the classical value of correction to the first 
moment of g\ is reproduced if we take into account the leading term of the polarization 
vector rj = p\ q /ra q . However contrary to the paper we would interpret the result with 
C g \ = 5/3 in eqs.(|92|) as physical one. In our calculation we took the proton (quark) 
polarization vector in a general form. Using of the exact representation of the vector [[Dj]: 



1 



V 



\J{kiPi q ) 2 ~ m 2 m 2 q 



klPlq r 



Piq ~ m q ki (93) 



leads to exactly the same result (with C g \ = 5/3) for the first moment of g\. The 
second term giving the additional contribution ~ 2/3 cannot be neglected within the 
approximation under consideration. It can be easy understood from pure calculation of 



QED and electroweak corrections to the lepton current in DIS process [31], where the 
lepton polarization vector looks like hadron one (||). The explicit expression of this non- 
vanishing for m — > contribution within electroweak correction to the lepton current is 
defined by the formula from previous section. The more detailed description of this 



contribution treatment can be found in [32]. 



The result with C g \ = 5/3 was obtained under similar assumptions in the report pSj], 
The additional correction was discussed both within OPE (operator product expansion) 
and within improved quark-parton model in paper |50[| . It was shown in this paper that 
there is no contradiction between this result and classical correction obtained for massless 
QCD, if we carefully take into account the finite quark mass effects for coefficient function 
and matrix element. As it was reviewed in recent paper 52] the renormalization of the 



axial-vector current for massless quark completely suppresses the additional contribution 
obtained within massive approach (see section 4 of the paper). 

We note that Burkhard-Cottingham sum rule is held within taking into account mass 
effects at the NLO level. The target mass correction of the next order (~ m 2 /Q 2 ) was 



analyzed and was shown negative in the paper [51]. However, in the case of conserved 
currents (see []53|]) the polarized structure function <?2 with the target mass correction 
obeys not only Burkhard-Cottingham sum rule but and Wandzura-Wilczek relation too. 

Estimating the diagrams with one-gluon radiation we have also the result for QED 
radiative correction to hadronic current. The calculation within the quark parton model 
shows that QED correction is not so large, however there are at least two arguments to 
be taken into account it. First, the DIS structure functions are defined for the one photon 
exchange approximation as an objects including only the strong interaction. It means 
that transferring from the cross section to the structure functions we neglect these QED 
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effects, so their contributions have to be included to systematic error of corresponding 
measurements. Second reason is that there exist some measurements where the QED 
correction is important. One of the example can be calculation of a s by comparing 
theoretical and experimental values for the Bjorken sum rule jn]. In this case a s is 
suggested to be extracted from QCD corrections to the Bjorken sum rule, which is a 
series over a s , and at least five terms of the expansion are known. Simple estimation 
shows that the QED contribution is comparable with already fourth (or even third) term 
of the expansion. 

There are several papers devoted to calculation of QED and electroweak corrections 
to the hadronic current |23-25,3l|,37,^2|. As usual methods similar ours are used for 
that. Positive moment here is the keeping quark mass non-zero. From the other side 
there are some effects which are not considered within the QED calculations: taking into 
account the confiment effect in integration of soft region, consideration of the quarks 
as non-free particles. Thus the methods (see for example |4j|) developed for careful 
treatment of the QCD effects should be applied for the analysis of photon emission from 
the hadronic current. However, probably the best way to take into account the photon 
radiative correction to the hadronic current is further generalization of OPE technique for 
including of the QED effects. 



3 Corrections to the lepton current in other pro- 
cesses 

3.1 Diffractive vector meson electroproduction and code 
DIFFRAD 

The measurement of the cross section of the exclusive vector meson electroproduction 

e(fci) +p(p) — ► e'(fc 2 ) + v(p h ) +p(p 2 ), (94) 

can provide information on the hadronic component of the photon and on nature of 
diffraction. During several years the diffractive production of the vector meson has been 
the subject of the muonproduction ||73|-f75fl and electroproduction ||76|-[7^| experiments. 
Data analysis of these experiments is considerably affected by the QED radiative effects. 
At practice RC to the processes of electroproduction are taken into account using codes 
originally developed for the inclusive case (see ref. [79|, for example). 



The purpose of this subsection is to present the electromagnetic correction to ex- 
perimentally observed cross sections that was calculate by Bardin-Shumeiko approach in 
ref. gj|. 

We consider RC to three and four dimensional cross sections a = da/dxdydtdcfth and 
a = da/dxdydt. They are related as 

a = [dcPhCi . (95) 



o 



The four differential Born cross section can be presented in the form 



CO 



An 2 xy 



{y 2 a T + 2(l-y--y z Y)(a L + a T )j, (96) 
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where x and y are usual scaling variable, ot and <jl are differential cross sections of the 
photoproduction, -y 2 = Q 2 jv 2 and v is the virtual photon energy. 

For the observed cross section of the vector meson electroproduction we obtain 

<?obs = °oe S * nf (1 + 5 V r + 5 vac ) + o F . (97) 

As earlier the correction 8 vac comes from the effects of vacuum polarization by leptons 
and hadrons. 

The contribution of the infrared finite part can be written in terms of POLRAD 2.0 
notation [13, 14]: 



' ■ J ■ (98) 



7"min ^3 



g 4 J Q 4 



where R = v/f, f = 1 + t — kp^/kp and 2Mr maXim j n = ± as usual 5j = 1 for 

j = 1 and <5j = otherwise. The quantities 0jj depend only on the kinematic invariants 
and the integration variables r and The explicit expressions for them can be found 
in |16|| . ^-functions are kinematical ones without any integrals for unpolarized case, and 
only one integration are left in the case of polarization part. 

The dependence on the photoproduction cross sections is included in Ti\ 

T\ = S x a T , T-i = 2x(<r T + cfl)- 

The quantities (Jt,l have to be calculated for Born kinematics, but is calculated in 
terms of so called true kinematics. It means that they have to be calculated for the tilde 
variables instead of standard Q 2 , W 2 and t which have being already defined by eq. pgj). 

The important point is the dependence of the results on the maximal inelasticity v m . 
The inelasticity is calculated in terms of the measured momenta, so it is possible to make 
a cut on the maximal value of this quantity. If this cut is not applied the maximal 
inelasticity is defined by kinematics only. Below we give the formulae for v m in terms of 
the kinematic invariants 



4Q 2 v m = \^JX q - ^t 2 + 4m 2 Q 2 J - (S x - 2Q 2 + t q ) 2 - 4M 2 Q 2 . (100) 

The FORTRAN code DIFFRAD created on the basis of the presented formulae calcu- 
lates the lowest order RC to the diffractive vector meson electroproduction. The higher 
order effects are approximated by the procedure of exponentiation. The formulae for the 
cross section are given in a covariant form, so the code can be run both for the fixed target 
experiments and for the experiments at collider. 

Below we give numerical results for RC factor 

/ d(f> h a obs 

• ,J " ! - " — (101) 
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Figure 8: RC factor under kinematic 
conditions of HERMES; = 7.9 
GeV, t = -0.11 GeV 2 . 
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Figure 9: Observable and Born cross 
sections as a function of t under kine- 
matic conditions of HERA; y/s=300 
GeV; W=70 GeV; Q 2 =4 GeV 2 . 



obtained within the kinematic regions of recent HERMES and HERA experiments for 
exclusive p(770) electroproduction. 

The dependence on the v m cut in the region of HERMES kinematics is analyzed in 
fig. 8. The usage of the cut changes RC factor for small Q 2 and does not influence on RC 
for larger values of Q 2 . In the case of the cut usage we have to define v m as a minimum 
of the value of the cut and v m given by the kinematic restrictions (100). For fixed values 
of —t the cut influence on v m and RC till certain value of Q 2 only. 

The Q 2 -dependence shown in this figure is typical for the situation when the inelasticity 
cut is not applied. If this cut is used then the rise of rj when Q 2 goes down would be 
suppressed. Notice that the t-dependence is rather important too. Figj9| illustrates this 
last property. r\ crosses unity for —t ~ 0.25-0.3 GeV 2 and rises with increasing \t\. 
The large positive correction in this case is a result of the large phase space for photon 
radiation. The cut applied to the inelasticity (or E — p z ) can again reduce the value of 
the RC factor for large \t\. As a consequence of the i-dependence of r], the observed slope 
parameter also receives large RC. 

Now as application of the obtained result let us consider the process 

e(h)+p(p) ^ e(k 2 )+p(p')+ p(p V ), p^Tr + (p+) + ir-(p_), (102) 

that can be viewed as an off-diagonal Compton scattering analytically continued in the 
virtuality of the photon 7* to the vector meson mass j*p — > Vp and gives access to the 
whole set of the corresponding helicity amplitudes. 

The process ( |102| ) is analyzed experimentally by means of spin-density matrix elements. 
When measured, they give an indication of vector meson internal constituents motion and 
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Figure 10: The dependence of Born (dashed line) and radiative corrected (solid line) spin-density 
matrix elements on Q 2 under the kinematic conditions of Hl/ZEUS experiments. 



its spin-angular structure. The angular distribution of unpolarized vector meson decay is 
parameterized by fifteen matrix elements rf^r^ . For a long time it was believed, that 
their behavior complies with the s-channel helicity conservation (SCHC) hypothesis, which 
means that the helicity of the virtual photon is conserved in the s-channel process "f*p — * 
pp. In this case ten matrix elements (which corresponds to the case when photon and 
vector meson have different helicities) are equal to zero. But in the recent measurements 
Tq has been observed to be non zero |80|-^1, what has been considered as an indication 
to SCHC violation. 

The procedure of the experimental data analysis is based on the correlation of the 
lepton scattering, vector meson production and decay planes, which are affected by RC. 
Hence it is topical to look at whether the measured Tq can, at least partly, be the 
result that RC coming from non-observed QED effects and real photon emission was 
underestimated. In any case in order to make the data processing of the corresponding 
experiments |8C- js!| to be consistent, RC should be taken into account. 



As it was shown in [18| the value of the QED corrections (Ar = r & s — rB orn ) to the 



matrix elements are defined by the following quantity: 



2tt 

/ d(f> h cos(p h 8(fa) 
o 

Jdfa (1 + S(fa)) 
o 



n = l,...,4, (103) 



where 6 (fa) = cf ohs jaQ - 1. 

For the majority of the matrix elements vanishing in the SCHC limit, radiative correc- 
tions turn out to be not greater then 1%. However there are two of them, Re r^o and rg , 
which RC appears to be substantial (see fig. |l0|). One can see, that corrections ARe r?o 
and ArQQ may reach ~ 20%. 
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The last result is interesting from point of view of the found SCHC violation: the 
radiative correction procedure reduces the observed effect. 



3.2 Elastic ep scattering and code MASCARAD 

Precise polarization measurements of nucleon form factors in electron scattering is an 
essential component of new- generation electron accelerators such as CEBAF [83]. This 
unprecedented precision requires knowledge of higher-order electromagnetic effects at a 
per-cent level. The purpose of this subsections is to analyze radiative corrections in elastic 
electron proton scattering and present computational techniques that could be used in 
experiments at Jefferson Lab and other electron accelerator laboratories. 
The observed cross section of the process 

e(k l ) + N(p)^e'{k 2 )+N( P2 ), (104) 

is described by one non-trivial variable, which is usually chosen to be squared of mo- 
mentum transferred. There are two ways to reconstruct the variables, when both lep- 
ton and nucleon final momenta are measured. In the first way it will be denoted as 
Qf = —{k\ — A^) 2 , and for second one it is Q\ = —{p 2 — p) 2 ■ It is clear that there is no 
difference between these definition at the Born level. However when the radiation of the 
lepton is only considered at the level of RC the symmetry is breaking. Here we will deal 
with both situations. 

In the first case the structure of bremsstrahlung cross section looks like 



da o f d 3 k 

~ ex. I — 

dQf J k 



J^lj2^ 2 (Ql)A (105) 



where /C is a kinematic coefficient calculatable exactly in the lowest order. It depends 
on photon variables. T 2 is a bilinear combination of nucleon formfactors dependent on 
Q\ only, which is function of photon momentum. Usually only the final momenta are 
measured in the part of the space. It is controlled by the function of acceptance .4, which 
is 1 or in depending whether the final particles make it in detectors or not. The integral 



(105) should not be analytically calculated for two reasons. The first one is dependence 



of formfactors on Q\. We avoid to use some specific model for them. The second one is 
acceptance usually very complicated function of kinematic variable dependent on photon 
momentum. 

For the second version of reconstruction of transfer momentum squared the structure 
of the cross section is 

da f d 3 k 

~ cr 



E^iQl) J (106) 



dQ\ - J k 

In this case the squared formfactor is not dependent on photon momentum and for 4tt 
kinematics (.4=1) this integral can be calculated analytically. In the experimental 
conditions at JLab |83]| , both final electron and proton were detected in order to reduce 
background. However elastic scattering kinematics was restored by the final proton kine- 
matics, while electron momentum was integrated over. Therefore, formalism related with 
Eq.( |106| ) applies for this case. 

Polarization effects are described by polarization four-vectors of the lepton (£) and 
proton (rj) which can be expanded over the measured momenta k\, p\ and p 2 - The 
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polarization vector of the lepton has been defined earlier by eq.(|5|) as well as polarization 
vector of the proton can be expand over four momenta of scattering particles 



77 = 2a v ki + b v q + c v (pi+p 2 ), (107) 

where the expansion coefficients defined by the different polarization states of the target. 
As a result the Born cross section for the process (104) can be written in the form 

dag 27TQ 2 1 



<IQ- * 2 Q' , 



where 



9? = 1Q\ Q% = ~^(Q 2 a v + (25 - Q 2 )c n ), 

Q ^ (109) 

9 ° = " ~ M2Q2) ' *° = W {2S ~ g2)(a " " 2bri 



The structure functions J- 2 are bilinear combinations of the nucleon formfactors dependent 
on Q\ only 

Tx = 4r p M 2 G 2 M , F 3 = -2M 2 G E G M , 

f 2 = Ul 2 ^±^M, T, = -M 2 G m Ge - Gm . (U0) 

l + T p l + T p 

where r p = Q 2 /4M 2 . 

According to |l^] the cross section that takes into account radiative effects within the 
leptonic variables can be written as 

da obs = da e S ^ (1 + 5 V r + <W) + da F . (Ill) 

Here the corrections 5i n f and 5 vac come from radiation of soft photons and effects of 
vacuum polarization. The correction 5vr is an infrared-free sum of factorized parts of real 
and virtual photon radiation. The infrared- free contribution of bremsstrahlung process 
can be presented as integral over three variables: inelasticity v = A 2 — M 2 (A = p+k\ — A^), 
r = kq/kp and angle 4>k between planes (q, k) and (ki, k2) 



or n i 7=1 



Here i? = v/(l + r) and „4 is integrated over the <\> acceptance function. The quantities 
Oij are similar to ones used in eq. (|l2|). However there they are integrated over cp^. We 
keep this integration because of possible dependence of acceptance function on this angle. 



The explicit expressions of the ^-functions in flll2j ) are discussed in Appendix of |1S|. 

Since within the hadronic variables the integration over the phase space of the real 
photon can be performed analytically, the observable cross section for this situation has 
more simple form 

da obs = IW + 4 ^ + C + S LWi]Fi + da . (113) 
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Figure 11: Born (solid line) and observed (dashed line) asymmetries vs Q 2 . No kinematic cuts 
on inelasticity were used. E = 4 GeV. 

The explicit expressions for Of and 8 el can be found in |20| |. 

On the basis of formalism presented in this subsections the FORTRAN code MAS- 
CARAD developed. This code uses Monte-Carlo methods to calculate RC to the observ- 
able quantities in polarized ep scattering measurements. 

Here we consider RC to some obsevables within two polarization measurements: 

1. The initial proton is polarized and the final electron is measured to reconstruct Q 2 . 
As a result there are four experimental situations for asymmetry definition: the 
target is polarized along (perpendicular) to beam or q (q = p2 — p). Corresponding 
polarization 4- vectors are denoted as t]l (t/t) or rf L (ry^). 

2. Polarization and momentum of the final proton are measured. Two polarization 
states should be considered: the final proton is polarized along (n' L ) and perpendic- 
ular (r]' T ) to q. 

It is natural that for the fist type of measurement we have to calculate RC within the 
leptonic and for the second type within hadronic variables. 

Both the spin averaged and spin dependent parts of the cross section (a u and a p ) can 
be presented as 

<f s = (1 + SK' P + (114) 

Both the factorized correction 5 and unfactorized cross section coming from 
bremsstrahlung process contribute to cross section. 



34 



5 U (%) 




5 L,T (%> 



u , GeV 



1.5 



0.5 



-0.5 




1 - Q 2 =1 GeV 2 

2 - Q 2 =3 GeV 2 



0.5 1 



2.5 



3 3.5 



Figure 12: Radiative corrections to the unpolarized cross section (left plot) and polarization 
asymmetries (right plot) defined in (JW. Solid and dashed lines corresponds to longitudinal and 
transverse cases. S=8 GeV 2 . 



Absolute and relative corrections to asymmetry can be denned as (see ( 114[ )) 

(i + *K + g£ 



A A; 



A; 



A 



A 



Ao ■ 
■ Ao 



0"n 



S p - 5 U 



o 



(115) 
(116) 



A i0 1 + 5 + 5 U 

where index i runs over all considered cases: i = L,T,qL,qT; S U)P = cr^ / o-q ,p . Here the 
correction 5 is usually large because of contributions of leading logarithms. However it 
exactly canceled in numerator of expression for correction to asymmetry. That is a reason 
why the correction to cross section can be large while the correction to the asymmetry is 
relatively small. 

Born and observed asymmetries are presented in fig. [ll]. There were no cuts for the 
missing mass. As a result due to unfactorizing properties hard photon contributions gives 
different contributions to spin averaged and spin dependent part of cross sections. For the 
longitudinal asymmetries 5 P > 5 U and there are positive contributions to RC. Opposite 
situation with the transverse asymmetries. 

For the second type of measurement let us define the relative corrections to the ob- 
servables in the current experiments: 

-1 



5 U 



1, 



- obs 

Up 








n obs 
_u 


< 







obs 



T obs 



1 



(117) 



The first quantity 5 U is the relative correction to the unpolarized part of cross section. 
The 5l,t give contribution to polarization asymmetries measured by rotating the polar- 
ization states of the initial protons. The correction to the unpolarized part of cross section 
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is presented in fig. |l~2|a. Its behavior is quite typical. For the very hard inelasticity cut 
(u m <C Q 2 ) the positive contribution is suppressed due to real bremsstrahlung, so there 
is only negative loop correction contributing to cross section. Different ending values for 
the curves corresponds to different kinematically allowed regions. 

Conclusion 

Some radiative effects in lepton scattering on polarized and unpolarized nucleons and light 
nuclei have been briefly reviewed. The processes of inclusive, semi-inclusive, diffractive and 
elastic scattering are considered. The explicit expressions for QED, electroweak and QCD 
lowest-order correction within Bardin-Shumeiko approach have been presented. Besides 
the next order correction in LO has been estimated. 

Basing on the presented formulae some FORTRAN codes have been created which 
allow one to provide the procedure of RC of experimental data in the current and fu- 
ture experiments as well. More popular of them are POLRAD Jl4| and RADGEN [O]. 
The version 2.0 of FORTRAN code POLRAD provides complete RC procedure in inclu- 
sive and semi-inclusive experimental setup. Monte Carlo generator RADGEN creating 
from FORTRAN code POLRAD 2.0 provides event simulation with taking into account 
radiative effects in various experimental setup on inclusive DIS. Codes DIFFRAD and 
MASCARAD calculates RC in diffractive and elastic processes. 

The obtained results and FORTRAN codes constructed on their basis are widely used 
in experimental collaborations at CERN, DESY, SLAC and TJNAF. 
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